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Abstract 

In this paper we discuss SU(N) Chern-Simons theories at level k with both 
fermionic and bosonic vector matter. In particular we present an exact calcu- 
lation of the free energy of the M = 2 supersymmetric model (with one chiral 
field) for all values of the 't Hooft coupling in the large N limit. This is done by 
using a generalization of the standard Hubbard-Stratanovich method because 
the SUSY model contains higher order polynomial interactions. 
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1 Introduction 



Three dimensional conformal field theories are of great interest in the context of AdS / CFT 
and dS/CFT correspondence [|TJ, [[|. An important feature is the duality |§ between the 
singlet sector of a vector model and a higher spin gravity theory [f§-[(|]. A typical example 
is the correspondence between the singlet sector of O(N) (U(N)) free/critical vector model 
and Vasiliev theory defined on AdS 4 with suitable boundary conditions [[| [7], [|. Evidence 
for this conjecture was provided by computing the correlation functions of the higher spin 



gauge theory PHTT 



Subsequently vector matter was gauged along with the addition of a conformally in- 
variant Chern-Simons term []12| , jL3[]. This ensures that all operators and states of the 
theory are gauge singlets. The corresponding dual higher spin theory also violates parity. 
Supersymmetric extensions of parity violating Vasiliev type higher spin theories have been 



recently constructed in [14]. Related works are | 



There are three main results in |L2||, where SU(N) level k Chern-Simons theory with 
a single fermion was discussed in the large N limit: i) calculation of the thermal 2-point 
functions and consequently a formula for the free energy as a function of the 't Hooft 
coupling and the temperature, ii) identification of the operator spectrum at large N with 
an infinite set of single trace gauge invariant currents of integer spin s and dimension 
s + 1 + 0(1/ N), iii) determination of the equation of the divergence of these currents, 
which is schematically given by 

d-j = jjjj + j^jjj. (1.1) 
Conservation and partial conservation of higher spin currents have been used to derive very 



general results about current correlation functions in [|T7], [L8| . 

The derivation of the thermodynamic results rested on the fact that i) the field theory 
was defined on R 2 xS 1 where one could choose the 'light cone' gauge[| in the non-thermal 
directions and ii) on dimensional regularization working in 2 + 1 — e dimensions. The 
choice of the gauge enables a representation of the path-integral in terms of bilocal fields 
which at the large N saddle point is the 2-point correlation function that solves the large 
N Schwinger- Dyson equations. 



This paper is a natural continuation of [12] by including not only a fundamental fermion 



but also a fundamental scalar (Chern-Simons- scalar theories were considered in |T3||). One 
of our main motivations is to study the M = 2 supersymmetric (SUSY) theory with one 
chiral field in the fundamental representation, which arises for special values of the param- 
eters in this class of theories. In this case, the conformal symmetry is exactly preserved 

El- 

The new features in the present study are the double trace Yukawa coupling and the 
triple trace (sextic) coupling that appear in the Lagrangian. To deal with them we develop 
a general method for carrying out path integrals in the large N limit that goes beyond the 
standard Hubbard-Stratanovich method ||20|| , [f2~l|1 that applies to quartic interactions. This 



general method is in fact valid for any polynomial action involving local as well as bilocal 



lr rhis gauge is explained in the main text, (X2. 
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fields and uses two auxiliary conjugate fields of a functional Fourier transform. At the 
large N saddle point these conjugate fields have the interpretation of the 2-point function 
(at large N) and the corresponding self-energy. The saddle point equations can be shown 
to be the Schwinger- Dyson equations at large N. 

Our exact result for M = 2 SU(iV) level k Chern-Simons theory with a fundamental 
chiral multiplet, on R 2 at temperature T, in the 't Hooft limit is the following. The 
self-energies of the scalar and fermion are 

S B (p) = cT 2 , Y, F ( V ) = f Ps + i( J ^--A 1 + P+ , (1.2) 

V Ps J 

where 7+ = -75(71 + ry 2 ), P+ = 775 (Pi - Wi), Ps = y/Pi+Pl (Pi>P2 are the components of 
momentum in R 2 ) and 



/ = ^ (log(2 cosh( + C )) - log(2 sinh(^))) , yf c = 2|A| ln(coth £). (1. 

The free energy density is 



3) 



67T 



yj\[ + Q <%log(coth|)J. (1.4) 



Note that the thermal mass c remains finite for all finite values of A and only diverges as 
c ~ In |A| as A — > 00. Similarly, the free energy is also finite and non- vanishing for all 
finite values of A. Its magnitude decreases monotonically from the free field value as |A| 
is increased and vanishes when |A| — > 00. These results are consistent with the conformal 



theory existing for all values of A. This behavior is also different from that in |jT2j| , where 
it was shown that for a fundamental fermion coupled to the CS gauge field the free energy 
vanishes at |A| = 1. 

The rest of this paper is organized as follows. In Section |2j, we present the exact solution 
of SU(iV) level k Chern-Simon theory coupling to fundamental scalar and fermion matter 
in the 't Hooft limit on R 3 . In particular, in Section |2.2.2| we present the general method 
for doing the large N path integral in the presence of higher order interactions using a 
functional Fourier transform. In Section Bl we solve the same system at finite temperature 
by analyzing it on R 2 x S 1 . Most of the analysis in these sections is carried out for a general 
theory consisting of one fermion and one scalar in the fundamental representation of the 
gauge group with varying values of the 't Hooft coupling, the Yukawa and sextic couplings. 
In Section |3.4| we specialise to the M = 2 case and present the results explicitly for that 
theory. Section £| is devoted to discussion and future directions. Appendices contain the 
conventions, technical details and consistency checks of the results in the main text. 



Note added: While preparing this manuscript we became aware of the work of O.Aharony, 
G.Gur-Ari and R.Yacoby P^j. We thank the authors for informing us in advance. We also 
mention that the thermal partition function involving only the CS plus scalar fields has 

[23. 



also been calculated in 
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2 Chern-Simons theory with fundamental matter in 
the 't Hooft limit 



2.1 Action 

In this section we present the exact solution of Chern-Simons theory coupling to fundamen- 
tal bosons and fermions with conformal symmetry in the 't Hooft limit. For this purpose, 
we specifically study SU(N) Chern-Simons gauge theory with level k coupled to a single 
fundamental boson and fermion ip on R 3 . The action is given by 

ie^^Tr(A„d„A p - %A P A V A P ) + D^Dy + ftfDrf 



S= d 3 x 

J I Ait " - — * - —' 

+ A 4 (^)(#) + A' 4 (^)(#) + A!' + (#)(#)) + A 6 (00) 3 



(2.1) 



Here Levi-Civita symbol e^ up is normalized as £123 = e 123 = 1. The covariant derivative 
acts on fields as follows (See Appendix [A| for relevant conventions). 

= - iA^Ty, = + %A a ^T\ (2.2) 

ZVA = - iA^, D$ = d$ + iA a ^T a . (2.3) 

The 't Hooft limit is defined by 

N, k 00 with A := f fixed. (2.4) 

Accordingly we define the coupling constants by 



A4 "^' A4 "2^' K ~Tk' Ae -(2^' (2 ' 5) 

where we have set k = j-. We keep the coefficients x 4 , x' 4 , x'l, x 6 of order one in the 't 
Hooft limit. Note in advance that the relevant coefficients to the self energy and the free 
energy are x^,Xq and the others do not matter in the 't Hooft limit. 



2.2 Exact effective action 
2.2.1 Gauge fixing 

The first important step to solve the gauge theory is gauge fixing. Since the gauge theory 
is being defined in R 2 x R 1 or R 2 x S 1 we fix the gauge degrees of freedom by using 
the following method. We first work with R 1 ' 1 x R 1 or R 1,1 x S 1 . R 1 ' 1 is a Lorentzian 
space-time. Here we fix the light-cone gauge 

A^ := -^={A 1 + A 2 ) = 0. (2.6) 
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Note that this gauge preserves Poincare symmetry in R 1,1 . In this gauge the Chern-Simons 
gauge interaction vanishes and the gauge field component A + = -75 (Al — A 2 ) is linear in 
the action. The action reduces to 

S = [ d 3 x \^-A a + d_A a 3 - 0(«9 2 + 2<9+<9_)0 + ^8^+ 



_2tt 

f A\J°L + A%Jl + A%A\$T a T h 



+ A 6 (#) 3 + A 4 (^)(#) + A' 4 (^)(#) + A: ((^)(^) + (#)(#)) 



(2.7) 



where J® = i(<j)T a d tl <j) — {d^^T 11 ^) — iip^^T "^. Here and hereafter we omit the total 
derivative terms by appropriate boundary conditions. Since A + is linear, we can integrate 
it out leading to the following constraint. 

^-d-A% + r = 0. (2.8) 
To solve this we Fourier transform to momentum space 

27T? 

A a 3 (p) = —J a _(p), (2.9) 

where p_ is understood using the principle value prescription p_ — > p_ + ie and 

J-(p) = -\ J ^((2q + P)-H-<l)T a <t>(q + p) + tH-qh-T a ij(q + p)). (2.10) 

Plugging this back we obtain an action with the gauge field integrated out, which is denoted 
by Si. At this point we analytically continue back from R 1,1 to R 2 . In terms of matter 
content, the action Si can be divided into three parts. 

Si = Sb + Sf + Sbf- (2-11) 
The first part Sb consists only of bosonic fields and is given by 

Sb = S ^M + ^~^) + N J ^^^^( P '^'*)^9i)x(-^<b) 

f d 3 P l d 3 P 2 d 3 qi d 3 q 2 d 3 q 3 . w w , , , 

+ N J WY 3 J2nf WY 3 (2^ j2^° 2{Pu P ' 9l ' fe g3)x(P ' gi)x(P ' q2U ~ Pl - P ' q3) - 

(2.12) 

Here x{Pi q) is a bilinear of the scalar fields with the gauge index contracted 

x(P,g) = ^0(^-g)0(^ + g), (2.13) 

where P represents the total momentum of the bilinear and q relative momentum. The 
coefficients Ci(P, q u q 2 ), C 2 (Pi, P 2 , q x , q 2 , q 3 ) are 

rm , 27rziV (-P + gi + q 2 ) 3 (P + + q 2 )_ 

C 1 {P,q 1 ,q 2 ) = — -, r , (2.14) 

k {qi - qst)- 

C 2 (Pi, P 2 , gi, g 2 , g 3 ) - ^ + - + - _ _ - + ^- + JV A 6 . 

(2.15) 
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The second part Sp is constituted only from the fermionic fields. The action in the leading 
order of large N limit is given by 

S ' = i / + /(§3(g3(g3^(- ? )7-^te + PW-)7sT"V.( 

(2.16) 

where we used ( |A.1|) and £i{P,q),£,-{P,q) are bilinear of the fermionic field with gauge 
and spinor indices contracted. 

tip, i) = ^(f - ?wf + g), UP, g) = ^(f - gh-i^f + g)- (2.17) 

The last term Sbf involves both the bosonic and fermionic fields. 

Sbf = N j 0- ^ 2 )3 2X 4 N X (P, qiM-P, to) + AS BF ( V , fj), (2.18) 
where ASbf(v>v) includes fermionic bilinear operators rj(P,q),fj(P,q) defined by 

V(P, q) = - q)^ + q), v(P, q) = - + q). (2.19) 

The precise formulas are available in Appendix [L).l| . For the large iV saddle points consid- 
ered in this paper, the classical values of these fields rj(P, q),fj(P, q) vanish. 

2.2.2 Method of functional Fourier transform for the large N limit 

In the previous subsection we have integrated out the gauge fields from the action ( |2.1| ) 
and obtained the action ( 2.11]) containing bosonic and fermionic bilocal fields. This leads 



to the complete elimination of the gauge field and the action is written entirely in terms 
of bilocal fields: 

N N N N 

^-(xo^fo), Y,^{xi)r(x2), j^ftwm EWso^rw 

m=l m=l m=l m=l 

(2.20) 

which are single trace operators, and kinetic parts like ip'y^d^ip, d^^d^, and some local 
terms. The appearance of bilocal fermionic and bosonic operators is a reflection of the fact 
that this theory contains higher spin fields. 

We now develop a general method for carrying out path integrals in the large N limit 
for arbitrary polynomial interaction.^ It is best to illustrate the main idea of the method 



2 Our method is an adaptation of a method developed for multi-trace unitary matrix model, see Section 
5 of p4 . 
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in a notationally simple example of bosonic fields. 

Spa-. — I d?xd^(j)d^(j) + J d i Xid 3 X2d 3 X3d 3 X4 : K 4: (xi 1 X2, X3, X4)((j)(xi)(p(x2))((j)(x3)(p(x4 i )) 

d 3 x 1 d 3 x 2 ■ ■ ■ d 3 x e K 6 (x 1 ,X2, ■ ■ ■ ,x 6 )((j)(xi)(j)(x2)){<j)(x 3 )(j)(x4))((j)(x 5 )(j)(x 6 )) + 



d 3 xd fl (pd> M <j) + V((jxj)). 
Introduce two conjugate fields a, fx such that 
1=1 VaS(a(xi, x 2/ 



—4{xx)<f>{x 2 )) 



VaVfx exp 



1 T , 



i I d xid x 2 fx(xi,x 2 )(a(x 1 ,x 2 ) ~ j^<P( x i)<P( x 2)) 



(2.21) 

(2.22) 
(2.23) 



We insert this identity into the partition function and integrate out the scalar fields <f> and 



V(pV(f)e~ 



VaVfx exp 



i / fx ■ a — V(Na) 



x / V<pV<f) exp 



% I d 3 x l d 3 x 2 d ll (j)d^(j) - fJ>(xi,x 2 )— 0(xi)0(x 2 )) 



VaVfi exp 



N ii Jx-a- V(a) - iVtr log(-<9 2 + ijx) 



where we set 



V(Na) = NV{a), fx{xi,x 2 ) = Nfx(x u x 2 ). 
After appropriate scaling by N the action is 

S ex = N Ur log(-<9 2 + ifl) + V(a) -i J Jx-a 

and the large N saddle point equations are 



5V 



5a(x,y) 
iot(x,y) 



+ ifx(y,x) = 
5 



5JI(x,y) 



trlog(-<9 2 + ift) = 0. 



(2.24) 



(2.25) 
(2.26) 

(2.27) 
(2.28) 



(2.29) 
(2.30) 



If V(a) is quadratic in a, the path integral with respect to a in Z reduces to Gaussian integral. 
This is the Hubbard-Stratanovich method and one is left with a Gaussian integral over /x. If V(a) is not 
quadratic, we have to work with both a and Jl. 
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From ( ggg ), ( ggg ) we obtain 



<*(x,y) = (x\ J M .~ \y). (2.31) 

— O z + Z/i 

We note that at the saddle point a is the propagator and iji = E is the self-energy to 
leading order in large JV. This equation can be written symbolically as 

(-d 2 + ijj)a = 1, (2.32) 

where the rhs is the identity operator with matrix elements (x\l\y) = 5(x — y). 
By using ( 2.29Q we find 

d 2 --^— )a = l. (2.33) 

These are precisely the Schwinger-Dyson equations after using large N factorization. 

J Vm^-y- [^y) m e- s ^] = 0. (2.34) 

This can be explicitly checked and it emphasizes the intimate relation between the existence 
of the saddle point equations and large N factorization. 

The free energy F ex = — log Z ex evaluated at the saddle point is 

F ex = N( [ tY\og(-d 2 +i(V)) + V{(a))-i J (2.35) 

where (a) , (Jt) are solutions of the saddle point equations. The above discussion could 
have been done as well in momentum space. 

It is important to note that the above discussion implicitly assumes the presence of an 
infrared cutoff and that the large N limit is taken before this cutoff is removed. In a sense 
N is the largest parameter in the theory. This method would fail if the infrared cutoff 
(volume of space) itself was of O(N). 



Let us apply this method to our case for which the action ( 2.11|) is given by 



s i= I 7^(c + <? 3 )<K-<?)0(<?) + * / 7^^/K-<?h A VK<?) + NV( x ,Zi,Z-,v,f)), 



(2.36) 



where 



V( X , £/, V,V) = J (2^3 (2^3 (^ ^(^ 9i, Q2)x(P, QiM-P, <fe) 



f d?P 1 d 3 P 2 d 3 q l d 3 q 2 d 3 q 3 
+ J (2^ WT 3 (2^ J2nf (2^ 2(Pl ' P2 ' fe gs)x(Pl ' giMP2 ' g2)x( " Pl 
f d 3 P d% d 3 q 2 8inN 

+ J j^j^j^ k ( qi - q2 )J- {p ^ l{ - p ^ 



The partition function can be written as 



V^VipV(j)V(j) exp 
d 3 q 



( 



q M + qM-q)Hq) 



(2tt) 3 

^(-q)YqMq) + nv( x , £/, ^)) 



(2.38) 



Now as illustrated above we introduce pairs of conjugate fields (a#, Hb),(®f,i, f^F,i), Hf,+) 
(ot v , n v ) (ctfj, such that 



J T>aBT>aF,iDaF-T>a rj T>a f 



~ 2/iF,+ • («F - 



Si a„ - - 



exp 



1 — 



N 



(2.39) 



where the factor 2 is inserted in front of \xf due to the fact that £/, £_ in ( |2.17 ) have the 
coefficient |. Inserting this identity into ( |2.38| ) we get 



i / jj, • a — NV(ctB, ocfj, ocf,-, oc^, a? 



x / VipVipV<f)V<f) e 



-S fr 



(2.40) 



where 



T>olD\l exp 



i J fj, ■ a = J VaBT>aFjVaF-Va v VafjViiB^f^Fj^f^F,-^f^ri^f^fj 



exp f z y (/i B • a B + 2/iF >+ • + 2/ip,/ • a F j + H v ■ a v + ' a v) S j 



and Sf ree is given by 

S free ~ 



(is + <?f )<K-<?)<K<?) + ^{-qh^q^iq) + i 



d 3 P 



(2.41) 



(-t*B(P,q)x(-P,q) 



(2tt) 3 V M V (2tt) 3 

2/i F , + (P, g)£_(-P, q) ~ 2n FJ {P, q) ~ ^(P, q)vi~P, <?) " ^{P, q)fj(-P, 



d 3 P d 3 q ,-.P . j.P 



(2vr) 3 (2vr) 3 v ^ v 2 



(4>{-- q )M^-q))Q{») 



0(f + q) 

m+q) 



(2.42) 



9 



which is the quadratic in elementary fields. Note that in ( [2.42 ) we have used 



(q 2 s + 4)K p ) + ¥b 
ip.fi 



i^q^5(P) + 7 + W,+ + ipF,i 



Hb = Nfi B , fi Fl = N/i Fl , n F+ = Nfj. F+ , fi v = Nfi v , fifj = Nfifj. 
By performing path integral for the elementary fields (j),(f),ip,ip we obtain^ 

VctV/j,exp N J p ■ ol — V(a B , a F j, &f-, ot v , a^) — STr log Q(fi 

Therefore we obtain the effective action as 

S e ff = N I STr \ogQ(fj,) + V(otB, Oi F j, a, F -, a v , ot^) — i p a 



(2.43) 
(2.44) 

(2.45) 
(2.46) 



In the rest of this paper we focus on calculations in the large N limit that are entirely 
determined by the saddle point values of the singlet fields. We will be interested in saddle 
point solutions that preserve both translational invariance and have fermion number zero. 
In this case the singlet fields become localf] 

(a B (P,q)) = (27r) 3 5 3 (P)a B (q), (a F> - s (P,q)) = (2n) 3 5 3 (P)a Frs (q), (a v (P,q)) = {a n (P,q)) = 0, 

(2.48) 

(ifJ. B (P,q)) = (2n) 3 5 3 (P)Z B (q), (i^(P,q)) = (27r) 3 5 3 (P)£ F , p (g), (i» v (P,q)) = (Z,- n (P,q)) = 0, 

(2.49) 

where 8 — 1,— and v — I, +. Using this it is simpler to derive the saddle point equations. 
The effective action Q2.46 ) takes a simpler form, as we explain below. 
The first term in ( |2.46| ) becomes 

STrlogQ(£) = STrlog( («2 + + E*(g)WP) 

(logfe 2 + ql + £*(<z)) - logdet(z 7 ^ M + 7 + £f,+ (<?) + £„,/(?))) 



(i^ qil + j+E F>+ (q) + E FJ (q))S(P) 



V 



(2vr) : 



(2.50) 



where we set V = (27r) 3 5 3 (0). The second term in ( |2.46| ) reduces to 
d 3 q 1 d 3 q 2 d 3 q 3 



V = V 

+ 



(2tt) 3 (2vr) 3 (2tt) 3 
d?qi d 3 q 2 87riA 



(2vr)3 (2vr)3 ( qi - q 2 ). 



C 2 (qi, 92, g3)«s(gi)ttB(?2)as(g3) 

-a F -(q 1 )a FJ (q 2 ) + [ -^-^-^-^8n\xia B (qi)a F j(q 2 ] 



(2vr) 3 (2tt) 3 



(2.51) 



4 Since ( 2.42 ) is quadratic in fields 4>, 4>, ip, ip, one can easily do the path integral by using Gaussian 
integrals. 



3 The corresponding formulas for the bilocal fields are 

( X (P,q)) = (2„) 3 S 3 (P) X (q), (f p (P,«)) = (2n) 3 5 3 (P)t s (q), (v(P,l)) = (v(P,l)) = 



(2.47) 



where 5 = 1,—. 
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where 



(gi + gs^fe + gs)- 



(2.52) 



£2(91, q 2 , qs) = C 2 (0, 0, gi, q 2 , q 3 ) = 4tt 2 A 2 I f- ;'" '-— + x 6 

In order to obtain ( [2.51 ) we have used the fact that the first term in ( |2.37 ) vanishes due to 
antisymmetry of the arguments inside the integral. Also note that the last term in (|2.37| ), 
^■S B p(a v , oifj) vanishes by using fl2.48|) . The third term in ( [2.46 ) which is written in ( p.41| ) 
becomes 

d 3 q 



fJL ■ OL — V 



(2tt) 3 

Collecting all the terms we obtain 
d 3 q 



-E B (q)a B (q) - 2E Fi+ (q)a F -(q) - 2T, FJ (q)a F j(q) . (2.53) 



S, 



eff 



NV 



(2tt) 3 
d 3 qi d 3 q 2 d 3 q 3 



+ 
+ 



(2tt) 3 (2tt) 3 (2tt) 3 
d 3 qi d 3 q 2 8-KiX 



(log(g s 2 + ql + Sfl(g)) - logdet(i 7 % + M?)) 
C 2 (qi, q 2 , q3)a B {qi)a B (q 2 )a B (q3) 



(2vr) 3 (2vr) 3 ( Ql - q 2 
d 3 q 

W) 



-a F -{qi)a, F j(q 2 ) + J 



d 3 q x d 3 q 2 



2tt) 3 (2tt) 3 

T> B {q)a B {q) ~ 2T, Ft+ (q)a F _(q) - 2E FJ (q)a FiI (q) 



8TrXx 4 a B (qi)a F j(q 2 ) 
(2.54) 



where 



Hp — ^f,+7 + + + Si? j3 7 3 + S^/J, with T lF _ = T lF3 = 0. (2.55) 

As we noted before (see Q2.31| )), at the saddle point of ( |2.54|) the singlet a and E 
fields have a simple interpretation as the propagator and the self energy, respectively. This 
follows immediately from the saddle point equations obtained by varying the action ( |2.54 ) 
with respect to the fields E B , Hp respectively.^ We obtain 

1 1 

a BW = o , ^ ,0 ot F {q) = - t— — 7-r , (2.56) 



q 2 + Z B (q)' 



where we have set 



a F = «i? i+ 7 + + cap -J + «f,37 3 + otpjl, with a F+ = a F ^ = 0. 



(2.57) 



-tr- 



Using (|2.56|) one can rewrite the remaining saddle point equations for a B , a F j and ctp,- 
as integral equations for £ fields in the following way. 

v / n n x f d 3 q 1 1 

S fl (p) 



-2ni\ 
2m X 



(27r) 3 (p-g)_ i^ qix + E F (q) 
d 3 q 1 7" 



-tr- 



(27r) 3 (p-g)_ ^ 7 ^ + S F (g) 



d 3 q 



47rAx4 



(2tt) 3 (2tt) 3 
d 3 q 



C2(p,q,q') + C 2 (q,p,q') + C 2 (q,q',p) 
1 



(2tt)V + ££(<?)' 
1 



tr 



(2tt) 3 2 7 ^ + E F (g)' 



g 2 + S B (g) g' 2 + £ B (g') 
(2.58) 



The equation 2(T,p t+ ap_- + Hpjoipj) = tr(Ef ctj?) will be useful for deriving the equation. 
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We refer to these equations as gap equations. 

As a check on our algebra we will rederive the gap equations (|2.58|) in two different ways. 
One is by using large N factorization on Schwinger Dyson equations of the fundamental 
fields in Appendix [FJ and the other by summing the planar diagrams contributing to 1PI 
self-energy diagrams, which will be studied in Appendix ||. 

By using the equations of motion (|2.56|) and ( |2.58| ) we can rewrite ( 2.54j) in terms of 
Eg, £_f as 



S, 



NV 
2 



d 3 q 



(2tt) 3 



M<? s 2 + <?! + Sfl(g)) - logdet(i 7 % + E F (q)) 

1 



3g? + 

47rAx4 

6 



1 

SbM + 2 tr 
d 3 g 1 
(27r) 3 g 2 



Sb(?) 



27^% + S F 
d 3 p 



(?) 



(27T) 



rtr 



+ E F (p) 



(2.59) 



This is our answer of the exact effective action in the 't Hooft limit. This means the 
effective action is exact in all orders of A. In a diagrammatic point of view, this effective 
action should encode the connected vacuum graph consisting of planar diagrams. We will 
check this point in Appendix [|]. 



2.3 Exact solution of the gap equation 

In this subsection we solve the gap equation ( [2.58|) . We note that the right-hand side of 
( [2.58] ) does not depend on p 3 . This means that Es(p), E^(p) are functions of momenta 



p + and p_. We first concentrate on the bosonic case, which is given by the last equation 
in ( |2.58D . Due to rotational invariance the self energy for boson only depends on the 



combination p s = ^j2p~^pZ. One can reach the same conclusion by noting that the bosonic 
self-energy E#(p) is real. This can be confirmed from the equation 

(0(-g)0b)> = (2tt) W(-g + p)a{p) = ^g^l^ (2 . 60 ) 

which follows from ( p.56|) . NowQ let us take the derivative with respect to p + on both sides 
of the last equation in (|2.58|). By using 



2tt5 2 (p) , -^—p- = and -^—p s = — , (2.61) 



dp + p~ dp + dp + p s 

we find that the first term on the right-hand side vanishes. The second term proportional to 
A4 is independent of the external momentum and therefore obviously does not contribute. 
Therefore we obtain 

J-^b(p) = 0. (2.62) 
dp+ 



7 This method for solving the gap equation was first explained to us in the context of the purely fermionic 
theory by S. Giombi. We thank him for discussions. 
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As a result £b(p) cannot depend on the momentum. Since we do not have any other 
dimensionful parameter, we conclude that 

S B (p) = 0. (2.63) 



Now we turn to solving the gap equation for fermionic case. For this purpose, it is 



helpful to compare with the purely fermionic case studied in |T^]. The difference between 
the purely fermionic case and our case is only the term with A4, which provides a constant 
shift independent of the external momentum in the equation for T^pj. £_f,+ remains the 
same as that for the purely fermionic case. The term proportional to A4 evaluates to zero 
by dimensional regularization 

where we have used = 0. Therefore the situation completely reduces to the purely 
fermionic case for which the solution was already obtained in [[L2|] to be 



= hp 8 I + WoP+7 + , fo = A, g = -A 2 . (2.65) 
We explicitly check that (|2.63|) , ( |2.65| ) satisfy the equation of motion (|2.58|) in Appendix 



3 Finite temperature 
3.1 Set up 

In this section we study SU(N) level k Chern-Simons theory with a fundamental matter at 
nonzero temperature, T. To this end we study the system on R 2 x S 1 with circumference 
f3 = 1/T. The boundary conditions are chosen to be periodic for the scalar and anti- 
periodic for the fermion. 

The important point is that the results for our computation on R 3 in the previous 
section can be easily adapted to the R 2 x S 1 case by replacing the integration over the 3rd 
component of the momentum with a discrete sum determined by boundary condition. For 
terms which involved internal scalar propagators we do replacement 




and for the part coming from an internal fermionic propagator, 

J V ' r P3-F ' nGZ y 

Note that the volume of space V is now given by V = V2P, where V2 is the volume of 
2-plane. 
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3.2 Exact self energy 



In this subsection, we solve the gap equation at nonzero temperature. Using the replace- 
ment rules described above the gap equations become: 



S F,/(p) 

S fl (p) 



-2-rri, 



qz-F 



d 2 q 



-tr- 



(2n) 2 (p-q)_ i^ + Z F (qy 
d 2 q 1 



7 1 \ - 

3f^J (27r)a(p-g)_"* 7 ^ + E f .(g) + 4 ^ 



-tr- 



d 2 g 



1 



«2 



47rAx4 



d 2 g 



(27T) 



£ 

d 2 g 



d 2 g' 



C 2 (p, g, g') + C 2 (q,p, q') + C 2 (g, g',p) 



(27r)V + £ B (g)' 

1 1 



■tr- 



(2tt) 2 z7^ + S F (g)' 



g 2 + S B (g) g' 2 
(3.3) 



Our strategy for solving these equations will be similar to that used in the zero temper- 
ature case in section [2.3| . We start with the third equation in ( p. 3D which determines 
The only dependence on the external momentum on the rhs of this equation is through the 
coefficients Ci{p, q,q') which are independent of the third component of the momentum. 
An argument similar to that in the zero temperature case shows that the dependence of 
the rhs on p+,p- also vanishes. It therefore follows that must only depend on the 
temperature and we can set 

E fl (p) = oT\ (3.4) 

where a is a constant with a non-trivial dependence on A which needs to be determined. 

Next we turn to the first two equations in ( |3.3|) which determine ^f,+- Once again 
the rhs of these two equations do not depend on p%. Rotational invariance in the 1-2 plane 
and dimensional analysis then allow us to write 



^f,i(p) = f(p)p s , Sf,+(p) = ig(p)p+, 

where f(p),g(p) are undetermined functions of p = ^, which is dimensionless. 



(3.5) 



3.2.1 Self energy of fermionic field 

Using the ansatz in ( |3.5|) the first two equations in (|3] 

d 2 q 1 f(q)q s 



take the form 



g{p)P+ = -47rA^ [ 

Q3-F 



(2 7 r) 2 (p-g)_(l + ^(g) + /(g) 2 )g 2 + g 2 



2 ' 



(3.6) 



Q- 



(27r) 2 (p-g). 



g(q) + /(g) 2 )?. 2 + 4 



+ 4vrAx 4 ^ I 

1 q 3 -B J 



d 2 q 



(27r) 2 g 2 + ctT 2 ' 
(3.7) 
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Now let us take the derivative with respect to p + on both sides. By using ( [2.6 1| ) we find 

£~ + mp+) = 5 jk (^ 9 (P) f ink)Pi^r (3 ' 8 

9 ,. M . . .l r 1 P_ 



= 4vrA^V-^-— — , P ,/~n 2 n 2 , 2 - ( 3 - 9 ) 



Combining these equations we find 





which implies 



{(f(p)Ps) 2 + g(M) =o, (3.10) 



/(p) 2 + ^(p) = 4, (3.11) 

where c is a dimensionless constant with a nontrivial dependence on A. From this relation, 
g can be obtained once we determine /. Therefore, we focus on solving the equation for /. 
Plugging ( p. 11|) back, we find 

(3.12) 



Let us first calculate the first term. First we do the angular part of the integration by 
using 

where <j) is the angular part of g_, g_ = q s e 1 ^, and 6{x) is a step function, 9(x) = 1 for 
x > 0, otherwise vanishes. Then we obtain 



Since this integral includes UV divergence, we regularize this integral as explained in 
Appendix 0. After the regularization, we find 

1 / — 2 I 

(1st term in (|3~lp ) = 4ttA x — - log (2 cosh( V:P )). (3.15) 

27rp 2 

Detailed calculation is available in Appendix p.3| . 

The second term in (|3.12|) can also be computed in a similar manner. Details of the 
calculation can be found in to Appendix |L).3| with the result 



(2nd term in ( gig )) = 4vrAx 4 x log(2 sinh(^)). (3.16) 
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The function \fa which appears in this formula is to be taken as the positive root. Also in 
obtaining this formula we are assuming a > 0. We will see below that this assumption is 
met by our solution in the M = 2 case for all values of A. Summing up these contributions, 
we find the solution of / and thus g as follows. 

/(p) = ^iog(2cosh(^±^))-x 4 log(2sinh(^))^). (3.17) 
9(P) = ^~f(p) 2 - (3-18) 
Here, again for convenience, we set p = f3p s . 



On determining the constant c in terms of a The solution for the self energy 
obtained above has two constants c and a . We now argue that the behavior of the integral 
equations (|3.6| ), (|3.7|) give rise to one relation between them. For this purpose, let us study 
the behavior of g(p) around p ~ 0. From (|3.17|) , ( |3.18|) , we find 



2A (log(2 cosh(^)) - x 4 log(2 sinh(^))) }' 



P 2 



+ 



(3.19) 



On the other hand ( |3.6| ), ( |3.17| ) imply that g(p) ~ Oil) around p ~ 0. To see this let us 
evaluate ( p.6|) by performing the summation over q% and the integral of the angular part 
of q. 



9(p)pa 



-4ttA 



-4vrA 



d 2 q 



tanh( v 2 



(27r) 2 (p-g)_ 2^q 2 s +cT 2 
p * q s dq s 1 tanh(^p^ 



f{q)q s 



Here we used ( |3TTT| ) , ( |TX23| ) and 



2tt 



2vr p_ 2y/q 2 s + cT 2 



Q{Ps ~ q s ) 



f(q)Qs 



(3.20) 
(3.21) 



2tt (p - g). 



P- 



(3.22) 



Therefore we obtain 



4ttA [*> q s dq s tanh( ^i +cT2 ) 



P 2 Jo 2tt ^2 + cT 2 



f{q)q s 



(3.23) 



Now ( |3.17|) implies that f(p) ~ i around p ~ 0, we find that behaves as g(p) ~ 0(1) 
around p ~ 0. Therefore, the singular term in terms of momentum in ( p.l9| ) has to vanish 
and we obtain 



(3.24) 



log(2cosh(^)) - x 4 log(2sinh(^)) 
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3.2.2 Self energy of scalar field 

Next we evaluate the self energy for the scalar field using the third equation of ( |3.3|) . Since 
this equation does not depend on the external momentum p, we can set p to zero. 



^ 2 = ^E/(^E/ ^ [g.(o, + cm, o, 4) + cm, rf, »)] g2 ^ ql2 ^ 

1 f d 2 q 1 
-4n\x 4 - > / 7— ^rtr ^^r- 3.25 

For convenience we denote the first term by Si and the second term by E2. First we focus 
on Ei, which takes the form, fl2.52Q , 

y 4 2 x 2 n+s \ 1 V f d2q d2q ' - - 

, 2 l r f d 2 q cPq> (q + q')- 1 1 



in 2 ' 



(27T) 2 (27T) 2 (q - q% [(^ + gj + E B (g)] [(^ + ^ + aT 2] " 

(3.26) 

Notice that the second term in Si vanishes due to the fact that the integrand is odd with 
respect to q, q'. Therefore, we need only to evaluate the first term, which can be obtained 
by applying the formula ( p.22|) twice, derived in Appendix p.3| . Thus we find 

2 



E x = 4tt 2 A 2 (1 + 3x 6 ) ( - -L l g(2 sinh(f))j . (3.27) 
Next we evaluate S2. 

E 2 = -4vrAx 4 i Y [ - . (3.28) 

Substituting T, F j(q) = f(q)q s into this we find 

S 2 = -47rAx 4 y (ii - x 4 log(2 s inh(^))/ ^ , (3.29) 

where we set 

_ 1 ^ f d 2 q log(2cosh(^|^)) 1 

1 P V J (2-) 2 ( 2 -^ly + ql + cT^ ° fi^J (2-) 2 C-^lf + ql + C T 2 ' 

(3.30) 

These integrals include divergence, so we have to regularize them to obtain finite results. 
We perform detailed calculation in Appendix p.3| (see ( p.24| ), ( p.22| )). The result is 

(log(2co S h(^))) 2 ! r c 

h = -~ 4^^< /o = -^>og(2co S h(^)). (3.31) 



17 



Therefore we find 



-l|^L(log(2cosh(^))) +2x 4 log(2sinh(^))log(2cosh(^))j. (3.32) 



Collecting these contributions, we find the solution for o as 



a — X 



ih(^)) J - 8x1 log(2 sinh(^)) log(2 cosh(^)) 



log(2cosh(^)) 



2 n 



(3.33) 



On the rhs the function y/a is to be taken as the positive root. 

Let us conclude this subsection by summarising the results obtained so far. The self 
energy for the boson is given in ( |3.4|) and for the fermion is given in ( |3.5|) with the functions 
f,g being given in ( [3.17D and Q3.18 ). These answers depend on two coefficients c, a which 
are obtained by solving the two equations Q3.24 ), ( |3.33| ). These expressions for the self 
energies are some of the main results of this paper. 



3.3 Exact free energy 

In this subsection we compute the free energy of the Chern-Simons-matter system on 
R 2 x S 1 . For this purpose, we divide the free energy density F into each contribution from 
scalar field, fermionic field and their interaction, denoted by F B , F F , F BF , respectively. 
They are given by 

S B (T) - S B (0) S f (T)-S f (0) S bf (T) - S BF (0) 
b = y , r F = — , -t bf = y , (3.34) 

where S B (T), S F (T), S BF (T) are respectively the contributions to the effective action com- 
ing from scalar field, fermionic field and their interaction on R 2 x S 1 

S F (T) = NV^J |^tr (- In + Mi)) + ^) iY% I Mg) ) , 0.36) 

~ 4it\x 4 NV 1 v f d?q 1 1 v f gp / 1 \ 

6 (3^J (2ny q i + Z B (q){3^ F J (2vr) 2 r Wp, + Sf(p) J ' 



(3.37) 



and S B (0), S F (0), S BF (0) are on R 3 . 

In what follows we compute F B , F F , F BF separately. 
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3.3.1 Contribution from scalar field 

First we compute Fb, which is given by 



B 



N 



-y 



d 2 q 
(2ny 



In 



q 2 + aT' 



aT 2 



3 (^f ) 2 + q 2 + aT 2 



d 3 q 



(2tt) 3 
(3.38) 

Here we already used the solution = aT 2 on R 2 x S 1 and = on R 3 . We compute 
this in the following way. 



In (q 2 ) 



N 



R 



d 2 q 
(2nf 



In 



2irn N 



+ 



d 3 q , ( q 2 + vT 2 . 



-y 

R L^i 



d 2 q 2 



(2tt) 3 
aT 2 



(27r) 2 3 \ (^f) 2 + g 2 + aT 2 



. (3.39) 



Since the first line in the bracket is the same as Casimir energy of free complex scalar 
fields with mass squared <tT 2 , it can be calculated as 



2 x 



/ W? log(l-e-^>) J^dy y log (1-e-y). 



(3.40) 



As first term in the second line of (|3.39| ) is linearly divergent, one can regularize it as 
explained in Appendix 



(2tt) 3 

We evaluate this as follows. 

d D q , f q 2 + oT 2 
^d 1o §(- 



d\ , q 2 + <yT 2 ^ 

log ( - ) —> hm 



r 



o 7 (2tt) 



lQ g( 55 i —o )• 



g" 1 + ?7r 



(3.41) 



(2*) 



g 2 + m 2 



d^q 



(2 



7T 



r (vf+^-vf+ 



(4tt)^ r(-i; 



(aT 2 



(m 2 ) : 



_ — (TVl - m 3 ) + 0(e) 

07T 



(3.42) 



where D = 3 — e is negative to make this integral convergent. Analytically continuing 
from some positive number e to zero, we evaluate this as — g^(T 3 <j2 — m 3 ), which becomes 
— g^(T 3 cri) under m — > 0. 

By using the formula ( p.21| ) one can evaluate the second term in the second line of 
( p.39j) , which gives 



3tt' 



a In ( 2 sinh 



(3.43) 



Collecting these terms we obtain 
NT 3 



^-o-i +2alog(2sinh(^)) + 6 J dy yln (l - e~")^ . (3.44) 

Note that the expression (|3.44|) is the same as the free energy density of Chern-Simons- 
scalar theory. 
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3.3.2 Contribution from fermionic field 



Next we concentrate on the contribution to free energy density coming from the fermions, 
F F . 



-^/(0 tr (- ln (^+4 o) (?)) + ^4 o: 



(3.45) 



where £p, E^? are the self energy of fermions on R 2 x S 1 , R 3 , respectively. We compute 
this as follows. 



N 



d 2 q 



log 



2n(n + 



2'\2 , „2 



d A q 



d 3 q , (ql + q 2 s +cT^ 



log 



(2tt) 3 
(2vr)3 2 gf + g 2 j " 



«i + <? 2 



(2tt) 3 

d 2 g 1 E Fj/ (g) 2 + cT 2 



log {ql + g 2 + cT 2 ) 



( 27 r)2 2 ( M^i) )2 + gs2 + cT2 



(3.46) 



The first-line in the bracket describes Casimir energy of free spin half fermions with 
mass squared cT 2 , and thus it can be evaluated as 



-2 x 



d 2 q 



log 1 + e 



-(g 2 +cT 2 ) 



I! 

7T 



/■oo 

/ ^^/logfl + e^). (3.47) 



The first term in the second line of (|3.46|) was already evaluated in the scalar case. The 

3 3 

result is -^r^-. 

By substituting the solution of Hp, the second term in the second-line of Q3.46| ) can be 
written as 



1 /2A 



J 2 -2x 4 log(2sinh(^-))/ 1 



2 V/3 

where I a (a = 0, 1, 2) are given by 



x 4 log(2sinh(^-)) ) + cT 2 }I 



(3.48) 



d 2 q 



log f 2 cosh(^p^; 



(27T) 



- 27r ( n +2^ 2 i „2 



- +g 2 + cT2 

These are calculated in Appendix p.3| , ( p.24|) . The result is 

(log (^2cosh(f)j ; a+1 

Ja ~ 2vr/3(a + l) 

20 



(3.49) 



(3.50) 



The final term in ( ft .46 ) is given by 

A 2 f d 3 q q 2 



2 J (2vr)3gf + ^ 



(3.51) 



which can be computed in a similar manner to ( p,12|) . One finds that this term vanishes 
as ( gXTg) does. 

Collecting these contributions we obtain 



y/c -S Qdyy\og{\ + e~ y ) - 2A 2 (log(2 cosh(^))) - ^clog(2 cosh(^)) 



NT 3 

+ 6x 4 A 2 log(2sinh(^)) (log(2 cosh(^)) J - 6A 2 x 2 (fog(2 sinh(^))J log(2 cosh(^)) 

(3.52) 

Note that when x 4 = Fp reduces to the free energy density of the purely fermionic 



system studied in [|T2j] . 



3.3.3 Contribution from interaction piece 

Finally we compute the contribution coming from the interaction between scalars and 
fermions. 



BF 



4ir\x 4 N ( 1 



d 2 q 



-y 



(2vr) 2 g 2 + g 2 + aT2/3^7 (2tt)^ \i^ Pfl + Z F (p) 



d 2 p 



tr 



d 3 q 



d 3 p 



tr 



(2tt) 3 ql + q 2 J (2vr)3 ^p^ + ^^p) 
We first compute the first-line. By using the formula ( |D.21| ), we obtain 

4vrAx 4 A r ( 1 , . , , v'^.A. . 1 \^ f d 2 p 2£ F ,j(p) 



(3.53) 



6 

N ( 1 



^log(2smh(f)))xI S / 



(27T) 2 ( M^i) )2+ ^ + cT2 



„log(2sinh(-^— )) x S 2 , 



(3.54) 



where we used ( |3.28| ) to obtain the second line. S 2 can be evaluated in a similar way as in 
fl3.32j ), so we get the first-line as 

N . . . 4A 2 x 4 / 



127T/3 



log(2sinh(— ))x- 



/3 2 



log(2 cosh(-^-)) ) - 2x 4 log(2 sinh(-^— )) log(2 cosh(-^-)) 

_ — 

(3.55) 



We note that the second line in (|3.53|) vanishes (see (p.4j , p.5l) for details). Thus we obtain 



NT 3 x 4 X 



BF 



^ ^log(2sinh(^)) ^log(2cosh(^))J -2x 4 (^log(2 sinh(^))^ log (2 cosh(^)) 

(3.56) 
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Summing up all terms and simplifying we obtain the free energy density as follows. 
NT " r -v^ 3 + v^ 3 + 2alog(2sinh(^)) - ^clog(2cosh(^)) 



6n 



- 2X 2 log(2 cosh(^-)) <^ log(2 cosh(^-)) - x 4 log(2 sinh(^— )) 



6 / dyy \og(l - e-y) - 6 / dyy \og(l + e~v) 



(3.57) 



The functions y/c and yfo which appear on the rhs stand for the positive values of the 
square root function. 

Th expression for the free energy density in (|3.57|) , which is valid to all orders in the 
't Hooft coupling A and the other couplings, x 4 and Xq, is one of the main results of this 
paper. Note that the coefficients c, a which appear in (|3.57|) can be obtained in terms of 
the couplings and the temperature from (|3.24|) and (|3.33|). 



3.4 Supersymmetric case 
3.4.1 TV = 2 case 

So far our analysis has been carried out for general values of the couplings £4, xq defined in 
(p.5|). One of our main motivations is to study the M = 2 SUSY theory in our setup and 
we turn to applying our general results to this theory now. For the M = 2 theory x 4 , x 6 
satisfy the relation^ 

£4 = xq = 1 (3.58) 

in the large N limit. More details on the M = 2 theory can be found in appendix F.l| . 

The self energy of the boson and the fermion was given in (|3.4|), (|3.5|), (|3.6|), (|3~7j), 
in terms of the constants a and c which are determined by ( |3.33|) , ( |3.24|) . The equation 
( |3.33| ), which determines a, for the choice of couplings ( |3.58| ) reduces to 



a = 4A 2 ( log(2 sinh(^-)) - log(2 cosh(^)) ) . (3.59) 



On the other hand (|3.24|) becomes 



c = 4A 2 ( log(2 cosh(^)) - log(2 sinh(^)) ) . (3.60) 



From these we findP] 



'c 



y/^= ^f c = 2|A|log(coth^-). (3.61) 



8 The J\f = 2 theory has two additional couplings x' Al x'l defined in (2.5) which take values x' A = 1, x'l 
in the large N limit. Their effects are suppressed in the large N limit, however. 
9 \fo, yjc in the solution stand for the positive root. 
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(a) 



(b) 



Figure 1: In Fig. [L(a) , \fc and \/a are shown as a function of A by undotted and dotted 
lines, respectively, when M = 2. They are degenerate and finite in all region of A. In 
Fig. |l(b)| , the free energy density normalized by —NT 3 is also finite in all region of A. 



This implies that the thermal masses of the scalar and fermion become the same in M = 2 
supersymmetric case. Correspondingly there is also a simplification in the free energy 
density which in this case becomes^ 



NT 3 
6tt 

NT 3 
6tt 



|A| 



6 / dyy log(coth -) 



> 



IAI 



6 



dyy log(coth- 



(3.62) 
(3.63) 



As seen from Figure |l(a)| , |l(b)| , the parameters a/c, y/a and the free energy density are 
all finite in the whole region of |A| G [0, oo). This is to be contrasted with the case of a 
fermion coupled to the CS gauge field [12) where the thermal mass diverges and the free 
energy density vanishes at |A| = 1 We will comment on this difference more in Section |3.5j 

The asymptotic behaviors as |A| — > 0, oo are also worth commenting on. Around |A| ~ 
we see from (|3.61|) and ( |3.62|) that a/c and the free energy behaves as 



a/c ~ — 2A log A(— log A) lo s A 
F - 



NT 3 /21 3 

tC(3) + 8A 2 (-logA)^(logA) 2 (31og2(l + logA) - log A) 



6tt V 2 

Around |A| ~ oo, \fc diverges and the free energy density vanishes as 



(3.64) 
(3.65) 



a/c ~ log 4A — log log 4A + 
NT 3 (log(4A)) 3 



log log 4 A 
log4A ' 



6/T 



A 



(3.66) 
(3.67) 



10 Note that we obtained the same results for c ( 3.61 ) and the free energy density ( 3.62) b y momentum 
cut-off r egula rization. We checked that the divergences appearing in the gap equations ( 3.61 ) and the free 



energy ( 3.62| ) in this regularization scheme cancel. See also |25) for cancellation of divergence in Af = 2 



case. 
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(a) 



(b) 



Figure 2: In Fig. |2(a) , we plot the yfc and yfa as a function of A by undotted and dotted 
lines, respectively, when w = 0. They are finite when < A < 2 and divergent at A = 2. In 
Fig. |2(b)| , the free energy density divided by —NT 3 is shown as a function of A with w = 0. 
—F > when < A < 2 and F = at A = 2. 



3.4.2 JV = 1 



case 



Another interesting case to consider is that of M 
realized when the coefficients (|2.5|) are given by 



1 SUSY0. M = 1 supersymmetry is 



w 



£4 



X A 



w, X A 



W — 1, Xq 



w 



(3.68) 



where w is a constant, which appears as a coefficient of superpotential ( |F.36| ). More details 
on the Chern-Simons-matter action of this theory can be found in Appendix F.2| . The large 
N limit of this theory is therefore characterized by two parameters, A and w. 
Let us now briefly consider the behavior for varying values of w G [0, 1]. 



w = In this case the M = 1 superpotential (|F.36|) vanishes. The graphs of a/c, a/cx are 
shown in Fig. 2(a) and that of free energy density (with some normalization) is in Fig. 2(b) . 

From Fig. p(a)| , one can read off that yfc and y/a monotonically increase up to A = 2 
and are divergent at that point. Since y/c and \fa correspond to the thermal masses of the 
fermion and scalar respectively we see that these masses diverge at A = 2. Correspondingly 



the free energy density with minus sign, which is shown in Fig. [2(b)| , is a monotonically 
decreasing function of A and vanishes at A = 2. This behavior of the system is consistent 
with the following picture. As the coupling constant A increases in value the thermal 
masses grow and the finite temperature correlations become shorter ranged. Finally, at 
A = 2 the thermal masses diverge causing the effective degrees of freedom to vanish and 
thus leading to a vanishing free energy density. 



< w < 1 In this region, a/c, y/a and the free energy behave qualitatively in the same 
manner as the w = case. Both a/c and \fa monotonically increase with respect to A and 



blow up at some A 



A cr , and the free energy, with a minus sign, decreases monotonically 



11 Application to M = 1 is first suggested by S. Kim. We thank him for discussion. Issues related to 
conformal invariance in this case will be briefly discussed in Section H 
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(a) 



(b) 



Figure 3: In Fig. |3(a)| , yfc and yfa are plotted as a function of A by undotted and dotted 
lines, respectively, when w = 0.5. They are finite when < A < 4 and divergent at A = 4. 
In Fig. |3(b)| , the free energy density divided by —NT 3 is shown. — F > when < A < 4 
and F = at A = 4. 



in term of A and vanishes at A = A cr . We show representative graphs of these functions in 
Fig. |3(a)| , |3(b)| when w — \ and A cr = 4. 



The w — 1 case corresponds to the J\f = 2 theory. As w approaches this value one find 
that A cr goes to larger values eventually becoming infinite when w = 1, in agreement with 
our discussion of the M = 2 theory above. 

3.5 Additional comments 

In this subsection we comment further on the results obtained above for the solution to the 
gap equations and the free energy at finite temperature. We saw in the previous subsection 
that for the M = 2 case the magnitude of the free energy \F\ monotonically decreases from 
its free value and vanishes in the limit when A — > oo. The thermal masses for the fermion 
and the boson are also well defined for all finite values of |A| and diverge at |A| — > oo. 
This behavior is different from what was found in the study of a single fundamental 



fermion coupled to the CS gauge field ||12|| . In that case the free energy vanishes at |A| = 1 



with the thermal mass also diverges at that value of |A|. This behavior suggests that 



the conformal theory stops existing beyond |A| = 1. An explanation was given in jT2 



as follows. First, consider a theory in the UV which contains a Yang Mills term in its 
Lagrangian besides also containing a CS term. The YM term dominates in the UV and 
this theory can be regulated using any standard regulator. Now starting with such a theory 
suppose we match in the IR to a theory which only contains the CS term for the gauge field. 
The YM coupling, g\ M has dimensions of mass so it would natural to do this sufficiently 
below the gfj V mass scale. On carrying out this matching one finds that the level of the 
CS term in the IR theory is shifted compared to its value in the UV theory. This shift was 



calculated in [26] and found to be 

\k IR \ = \k uv \+N. (3.69) 

In the UV theory, since the UV behavior is that of a conventional YM theory, kjjy can take 
all values and therefore \Xuv\ £ [0, oo). From ( |3.69|) it then follows that A = ^- satisfies 
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the relation 

|A| = - r 7;' < 1, (3.70) 









1 + 


^uv\ 



in agreement with the expectation coming from the free energy calculation that the con- 
formal theory does not exist for larger values of |A|. 

In fact an argument along these lines, now for the M = 2 theory leads to the conclusion 
that the conformal theory in this case should continue to exist for all values of A. In the 
M = 2 case if one starts in the UV with an additional Yang Mills term for the gauge 
field one must also add its supersymmetric completion which involves additional adjoint 
fermions. It turns out that when one matches to the IR theory without the SUSY Yang 
Mills term now there is no shift in the CS level. The gauge field results in the shift ( |3.69| ) 



but the two (real) fermions give rise to an exactly opposite shift making the net total shift 
in the CS level vanishP]. This conclusion is in agreement with our results for the M = 2 
theory. As mentioned above in this case the free energy does not vanish for any finite value 
of A (Figure l(b)| ) and the thermal mass of the fermion and boson also remain finite for all 



finite values of A (Figure |1 (a 



Another interesting case is one with M = 1 supersymmetry with the parameter w 
defined in ( |3.68| ) set to vanish. In this case the SUSY completion of the Yang Mills term 
requires us to add one real adjoint fermion. Integrating out this fermion results in a shift 
k k — N/2 so that the net shift, after including the effects of the gauge field, are 

N 

\km\ = \kuv\ + y (3-71) 
and the resulting value of the 't Hooft coupling is 

|A| = |A ^ y| , < 2. (3.72) 



When Xjjv — > oo we see that A — > 2 which is exactly the value we found in Figure [2(a) 
2(b) where the free energy vanishes and the thermal masses blow up. 



These observations serve as consistency checks on the results we have obtained. 



4 Discussion 

In this paper we have studied SU(N) level k Chern-Simons gauge theory with matter 
consisting of fermions and bosons in the fundamental representation. We calculated the 
self-energy for the matter fields and also obtained the free energy at finite temperature, 
in the 't Hooft limit. Our leading large N results, which were obtained using a saddle 
point method, are valid to all orders in the 't Hooft coupling. Although our emphasis in 
the study has been on the M = 2 theory with one chiral superfield in the fundamental 
representation, much of the analysis is general and would apply to other conformal field 
theories in this class as well. Our results show that the free energy, F, of the M = 2 theory 

12 In addition the supersymmetric multiplet also contains a real scalar. However integrating it out does 
not lead to a shift in the CS level. 
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is well behaved for all values of the 't Hooft coupling. Starting from the free theory, at 
A = 0, the magnitude of the free energy, \F\, monotonically decreases and finally vanishes 
at strong coupling, |A| — > oo. In addition the thermal mass of the boson and fermion, in 
the M = 2 theory, also remain finite for all finite A. The existence of the solution for all 
values of A is consistent with the exact superconformal symmetry of this system for all 
values of N and k. We have discussed the relation of our system to UV theory in Section 



3.5 



Our analysis for the self energy and the free energy can also be applied to theories with 
M = 1 supersymmetry or no supersymmetry, as long as they are conformally invariant. We 
have not carried out a detailed analysis of when conformal invariance can be preserved in 



this larger class of theories. The results in [25 1 suggest that conformally invariant theories 



with M = 1 or no supersymmetry can indeed exist. We leave a more detailed analysis of 
this question for the future. 



Lorentz invariance Our calculations were based on the techniques developed in [12| and 



were carried out by choosing light cone gauge A_ = in the non-thermal directions and by 
using dimensional regularisation. The choice of gauge we have made is unconventional in 
Euclidean space and more study is needed to ensure that the resulting theory is consistent 
and Lorentz invariant. Some evidence was provided for this in |12[ especially in the case 



of the anomalous dimension of the operator ipip which was calculated up to 2-loops in 
perturbation theory in various gauge. Our calculations which explore a larger class of 
theories add to this evidence. It is worth noting that, in the theories we study, the poles 
in the zero temperature propagators occur at p 2 = 0, which is Lorentz invariant. 

Higher spin currents It is also interesting to study the structure of the higher spin 
currents in our theory. The higher spin currents for such theories with vector matter have 
been discussed in [|17], |18[, and it has been shown that the resulting current algebra is a 
very powerful tool to determining correlation functions in the theory. These currents are 
also important in making a closer connection with dual Vasiliev theories. 

It is easy to see that in the theory we are studying, with one fermion and one boson 
in the fundamental representation, and with the double trace Yukawa couplings and the 



triple trace (jf coupling, ( |2.1| ), there are two currents for every integer spin s > 1 which 
are approximately conserved with a conservation law^] schematically of the form ( |1 . 1| ) . 
These currents can be built as follows. We start with say the current of spin s in the free 
fermion case, then replace all ordinary derivatives with covariant derivatives. The effect 
of the covariant derivatives, which do not commute, and the Yukawa and 6 terms in the 
Lagrangian results in the current not being traceless in general. After adding appropriate 
multi-trace terms traceless can be restored. This traceless spin s current is not conserved, 
again because of the non-commutation of the covariant derivatives and the presence of the 
Yukawa and 6 couplings. The resulting terms which violate current conservation are of 
double trace and triple trace type and are suppressed by appropriate powers of N as shown 
in (|1 . 1|) . A similar construction can be carried out starting with the spin s current in the 



free boson case resulting in the two approximately conserved currents for each integer s > 1 



13 There is an exactly conserved current for spin s = 1,2. 
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that we mentioned above. In addition, this is important in the analysis of SUSY theories, 
there are partner currents which have half-integer spin. 

It would be very interesting to explore current equations (analogues of ( |1 . 1| ) ) for the 



ABJ type theories |T4], |27|]. For example in the case of M = N and the adjoint repre- 
sentation of U(N) the gauge invariant currents of higher spin s > 3 are classically not 
conserved. The divergence of the current is not suppressed by powers of 0(1/N), as can 
be easily deduced by large ./V power counting. This has the implication that all higher spin 
(s > 3) gauge symmetries of the bulk theory are explicitly broken presumably by boundary 
conditions. 



ABJ duality Adding more flavors and the exploring theories with more supersymmetry 
is interesting. Particularly interesting is the case with J\f = 6 SUSY [27]. It has been 



argued that the Af = 6 theory has the duality |28| between the theories with gauge groups 
U(M + l) k x U(M)_ k and U(M) k x U(M + k- Z)_ fc . For the case where M = 1,1 = N this 
reduces to a duality between a U(N + l) k x [/(].)_& theory and a t/(l)fc x U(l + k — A)„ fc 
theory. The 't Hooft coupling of the non-Abelian groups on the two sides of this duality 
in the large N limit are A = x- and A = 1 — A respectively. The free energy of the two 
theories at finite temperature should be equal from this duality, once we also account for 
the change N — >■ k — N. The U(l) theory is in the large flavor limit here. 

We speculate the following two reasons to explain why the free energy formula presented 
in this paper does not have the duality: 

i) The presence of an additional U(l) and many flavors in principle leads to modified 
gap equations. Even though the saddle point of the SU(N) theory continues to solve 
the new gap equations there are in principle other solutions, whose inclusion in the path 
integral may restore the duality. 

ii) The large N saddle point analysis presented in this paper may be rendered invalid 
due to infrared divergences in sub- leading orders in However this possibility seems 
remote because both fermions and bosons have 0(1) thermal masses. 

iii) The use of manifestly lorentz invariant gauge may also throw light on this issue. 

Some other interesting problems to pursue include comparing our calculation with the 
high temperature limit of the theory on S 2 x S 1 where a Gross- Witten-Wadia type phase 
transition [[29|, |3(| is expected at a (dimensionless temperature) T ~ A 1 / 2 [Q. It will also 



be interesting to calculate the partition function of the SU(N) theory on S . Other gauge 
groups are also interesting to consider, especially in the case with Sp(N) gauge symmetry 
in the context of dS/CFT in |2], [31], |32[ . This has been independently suggested by [[33]. 
Finally, it would be worth exploring connections with Vasiliev theories in more detail. 
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A Conventions 

In this appendix, we collect our convention used in this paper. 

Normalization on gauge group The SU(N) gauge field will be expanded by the 
generators of the gauge group T a as = A^Ta, where a runs from 1 to N 2 — 1. T a is 
normalized so that tr(T a T b ) = S ab . Note that 



E( T XCH« = « - hy*. (a.i; 



N 

a 

Spinor contraction We contract spinor indices for fermion bilinear in the following way. 

iiA X ■= ^A a 0X fi = e ai ^A a pX ^ (A.2) 
where A is arbitrary two by two matrix and e a p is e-invariant tensor. 

Gamma matrices (Euclidean) We choose the gamma matrices 7^ as Pauli matrices. 

1 A . „,2 2 _ / ~i\ . „3 3 /I 



tW=I J ; 7W= ; • ; 7 w= 3 ). (a. 3) 



Coordinates In the main text, we have considered Lorentzian space-time R 1,1 from (|2] 
to ( [2. 11|) . x ± , A T ,d T ,p T represent 



± x x ±x 2 A ± ±A 2 d 1 ±d 2 Pi±P2 fK 

x = p— , A T = -= — , d T = p T = 7 ^-. (A.4) 

y/2 V2 V2 V2 

After (|2.11|) , we consider Euclidean space R 2 and x ± , p T represent 

± x 1 ± ix 2 pi ± ip 2 . . 

We denote the absolute value of 2-plane momentum as p s . That is, 

: = Pi + = 2 P+P-- ( A - 6 ) 
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Fourier expansion For a complex field we can define by Fourier expansion. 

On R 3 , can be expanded as 

On R 2 x S 1 , can be expanded as 

$ (*) = ^E / ^ ipx ^i ( A - 8 ) 

where p 3 = p 3>n is determined by boundary condition of the field $ with respect to S 1 . 

Gaussian integration The Gaussian integration in momentum space is done as follows. 
For a complex scalar field and a positive function with respect to momentum a(p), 

V<t>V&- J ^* ( - ria(pW(p) = (Deta(p))" 1 = e~ V 1 (& loga(p) , (A.9) 

where we set V = (27r) 3 <5 3 (0). For a fermionic field and a positive-definite two by two 
matrix M(p), 

~V^e~ J ^ { - p)M(pmp) = DetM(p) = e V f (A.10) 



B Schwinger-Dyson equation 

In what follows we shall derive Schwinger-Dyson equations for scalar and fermion field and 
confirm that they reproduce the gap equations obtained as a saddle point of the exact 
effective action. We will follow p3 . 

Schwinger-Dyson equation for scalar fields First we derive Schwinger-Dyson equa- 
tion for the scalar field. The Schwinger-Dyson equation can be derived via 

= / Dmj^^ , (B.l) 



where S\ is given by (|2.11 ). Since we are interested in gauge singlet sector, we contract 
the gauge index. 

= !/i*ZJ#£j-^ ('"TO) ■ (B.2) 
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Here we divide both sides by the total partition function. From direct calculation, the 
right-hand side becomes 



r.h.s. = N5 3 (q + p) 



2 3 



(27T) 



(2tt / 
d 3 qi d 3 q 2 
(2vr) 3 (2vr) 3 



(Ci(2gi - 2p, gi, g 2 ) + Ci(-2?i + 2p, g 2 , qi)) x 



x(g + 2gi - p, 9 + — -)x(-2g a + p, g 2 ) 



2 3 



(27T) 



<i 3 P c/ 3 g! d 3 q 2 d 3 q 3 
(2vr) 3 (2tt) 3 (2tt) 3 (2tt) j 



C 2 (2g! - 2p, P, g 1; g 2 , g 3 ) + C 2 (P, 2 ?1 - 2p, q 2 , q u q 3 ) + C 2 (2p - 2g x - P, P, g 3 , g 2 , gi) 
— ^ -(- 2 — 

X(g + 2g! - p, )x(P, g 2 )x(-2gx + 2p - P, g 3 ) 



x 



- 2A 4 iV 2 



2 3 /" rf 3 gi d 3 g 2 



(2vr) 3 J (2vr) 3 (2vr) 3 



£i(-2g 2 + p, gi)x(g + 2g 2 - p, 



-g + 2g 2 -p, 



(B.3) 



In the large N limit the leading term is dominated by factorized correlators. 



r.h.s. = N 



(2tt) 3 
2 3 f d 3 q x d 3 q- 



(2vr) 3 J (2tt) 3 (2tt) 3 



(Ci(2gi - 2p, qi , q 2 ) + Ci(-2g x + 2p, g 2 , q x )) x 



x(q + 2gi -p, 



-g + 2g x -p, 



(x(-2gi +P,g 2 )) 



2 3 f d 3 P d 3 qi d 3 q 2 d 3 q 3 



(2tt) 3 J (2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 

C 2 (2 qi - 2p, P, g 1; g 2 , g 3 ) + C 2 (P, 2g x - 2p, g 2 , g 1; g 3 ) + C 2 (2p - 2g x - P, P, g 3 , g 2 , g x ) 
-g + 2p-p/ 



x 



X(q + 2gi - p, 



-) )(x(P,g 2 ))(x(-2g 1 + 2p-P,g 3 )} 



2A 4 iV 2 - 



2 3 f d 3 q x d 3 q 2 



(6(-2g 2 + P, gi)) ( x(s + 2g 2 - p 



-g + 2g 2 -p, 



(2vr) 3 J (2vr) 3 (2vr) 3 



From (p|), (£171) , (12361) , (x(P,g)> and (£(P,g)> are given by 
1 „ 1 



Of— ) 



(B.4) 



(x(^?)> 



q 2 + Z B (q) 



(2tt)^(P), (£(P,g)> 



(2tt) 3 5 3 (P). (B.5) 



i^q^ + S F (g) 
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Using the above ansatz we can proceed 



r.h.s. 



N5 3 (q + p) 
q 2 + Z B ( q )[ 
1 



Sb(?) 



d 3 q 2 d 3 q 3 
(2tt) 3 (2tt) 3 



C 2 (q, q 2 , 93) + C2O2, g, g 3 ) + C 2 (q 3 , q 2 , g) 



g| + S B (g 2 )g| + S B (g 3 ) 



A 4 iV 



tr( 



(2tt)3 ^ 7 Mg'+E F (gO 



a) 



(B.6) 



Note that the quadratic term with respect to x vanishes and C 2 (p, q, r) is given by ( |2.52| ). 
Since this was required to vanish we obtain Schwinger- Dyson equation as 

d 3 q 2 d 3 q 3 



Ms) 



(2tt) 3 (2tt) 3 L 
d 3 q' 



X 4 N 



C 2 (q, q 2 , qz) + C 2 (q 2 , q, q 3 ) + C 2 (q 3 , q 2 , q) 
1 



tr 



(2tt)3 Hj% + ^ f W 
This is the same as the gap equation of the scalar field ( |2.58| ). 



Jg 2 2 + £ B (g 2 )gf + £ B (g 3 ) 
(B.7) 



Schwinger-Dyson equation for fermionic field We repeat the same procedure for 
the fermionic field. The starting point is the following identity. 

= J D^Di,^—^ (e-^V)) , (B.8) 

or by taking contraction for gauge index 



(B.9) 



In purely fermionic case this calculation was already done in [|12|]. The difference between 



their case and our case is one term coming from Sbf- The equation (2.27) in [[TJ 
corrected by a constant. 

N(2ir) 3 5 3 (p' + p) = i Pti r(ip m (p)$ m (p')) 



is 



+ 



2m 



d 3 r d 3 q 1 



k J (2tt) 3 (27r)V 
2ni f d 3 r d 3 q 1 



+ X 4 N 



(2tt) 3 (27r) 3 g 
d 3 q 



(2-it) 3 ip'^ + 
In the large N limit the factorization leads to 



1 + (MP ~ qW{-r)l 3 ^ m {r + qW{p')) 
1 3 (MP - q)i> n (-rh + Mr + qW(p')} 
l^r-rA^P +p-q)m)- (B.10) 



(2n) 3 5 3 (p' + p) 



i^Pn 

1 2ixi f d 3 r d 3 q 1 



■1 + (MP - q)r{-r)){l^m{r + qW{p )) 



+ 



il^Pv k J (27r) 3 (27r) 3 g~> 

~1^ 2 -V I ff\s ft\ ^(MP ~ q)r{-r))l + (^ m {r + q)ffV)) 
ijVpn k J (2ti) 6 (27r) 3 q + 

~1^^ N [ ( f%-> ^ 7 o fr(p'+P-#))- (B.11) 
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Plugging the fermionic and bosonic propagator into the above gives the fermion self energy 
to be 



Si?(p) = — 27riA 
+ X A N 



1 



7 



d 3 q 
dP 1 



1 



1 



i^q^ + S F I (p-q) 



(27r) 3 p 2 + E s (p)' 
which is the same as the gap equation for the fermionic field ( 2.58|) . 



(B.12) 



C On regularization 

In this appendix we have a comment on the dimensional regularization we used in this 
paper. In our analysis, we encounter momentum integrals which includes UV divergence. 
They are of the form 

j^npM) or l^j^F^M), (CI) 

where F(p s , \p\) is a function of p s = \pp\ + p| and \p\ = \Jp\ + p\ +p§. To regularize UV 
divergence, what we do is to change the dimension from 3 to D = 3 — e in the integration 
and also momentum according to it. 

/(§^(P,|P|) « ^£/(g^(P,|p|), (C2) 

where |p| = ^/pf + p| + • • • + p 2 D - It is remarkable that we do not touch the part of 
2-momentum p s under this process. It is clear that the UV divergence is cured if e is 
greater that some positive number. Performing integration, we obtain some function of e. 
Analytically continuing from some positive number e to 0, we obtain what we wanted to 
evaluate. 

In case where the integral includes IR divergence, it can be cured by adding a virtual 
mass parameter in 2-plane by changing p s — > ^Jp 2 s + m 2 and setting m = after evaluating 
the integral. 



D Detailed calculation 

In this appendix we collect results which will be obtained after detailed calculation. In 
this appendix and the main text we will use the following integration formula without 
mentioning 

r d D q 1 

J (2n) D (<f + M 2 ) a (4vr)f T(a) M 2a ~ D ' 

which is available in a standard textbook of quantum field theory. 



T(a- 



D - 
2 ■ 



m i ^ 
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D.l ASbf 

AS BF (r],fj) is given by 



AS BF (ri,fj) = N 



f d 3 P d 3 gi d 3 q 2 ( 27iN(P + qi + q 2 ) 3 , 
J (2n) 3 (2n) 3 (2n) 3 { k (gi-fc)- l?K 



N\' 4 fj(P, qi)r)(-P, q 2 ) + iVAl'{r/(P, qi)r](-P, q 2 ) + r/(P, gi)r/(-P, q 2 )} 

Ar /• d 3 ^ d 3 p d 3 d 3 d 3 2n 2 N 2 
N / — - — , \ n , \„ , \ n — X 



(2tt) 3 (2tt)3 (2tt)3 (2tt)3 (2tt) 3 k 2 
tr (^P 1 ,q 1 h-v(P2,q2)v(-Pi - P 2 ,q 3 h- ~ v(Pi,qihMP2,q2h-V(-Pi ~ P2,q 3 )) 

+ i(P 2 - P 1 - 2 qi - 2q 2 - 4g 3 )_r/(Pi, qi)j-Ti(P 2 , <fe)x(-Pi " ^2, g 3 ) J • (D.2) 
D.2 Check of the solution of the gap equation 

In this appendix we confirm that the solution presented in Section |2.3| satisfies the gap 
equation ( 2.58Q . That is, the equation we should like to check is 

d 3 q d 3 q' 1 1 



3iV 2 A f 



+ 8tt 2 A 2 



+ 4tt 2 A 2 
- iVA 4 



(27r)3(27r)3g2 q >2 

d 3 q d 3 q' r(p + q')-{q + q') 



1 1 



(2tt)3 (2tt)3 l(p - q')_(q - q')_\ q 2 q' 2 

d 3 q d 3 q' \(p+ q)-(p + q')-l 1 1 

(27r) 3 (27r)4(p- g )_(p- g ')- 
d 3 q 



q 2 q' 2 



tr- 



(2tt)3 iqK lfl + Z F ( q )- 

We compute the integrals by dimensional regularization explained in Appendix 0. 
First we focus on the first term, which is given by 3N 2 X 6 J 2 with M = 0, where 



(D.3) 



Ji 



d D q 



{2n) D q 2 + M 2 ' 



By using ( p.l|) we can evaluate this as 



1 



(47T 



§ M 2 ~ D 



M 
An 



0(e) 



(D.4) 



(D.5) 



where D is greater than two for the integration to be convergent. J\ — >■ — ^- after analytic 
continuation with e — >■ 0. Since J\ vanishes under M — 0, the first term also does. 

Next we compute the second and third integrals. We consider the second integral, 
8n 2 X 2 J 2 , where 



d 2 q d 2 q' r (p + <f)-{q + q')- 



?3 



<r 3 



(2n) 2 (2n) 2 l(p-q'Uq-q')^J (2n)^ q 2 + q 2 J (2n)^ q' 2 + q 2 ' 
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(D.6) 



Here we separate out the two light cone directions and remaining 1 — e directions. Let us 
perform the integration in 1 — e dimensions. 



J, 



A 



d 2 q d 2 q' r (p + q')-{q + q') 



(2tt)2 (2tt) 2 L(p - q'Uq - q%l \q s \^ \q' s \^ 



where we set 



A :-- 



r(i-i=£) 



0(e). 



To carry out the remaining integral, we use the following formula 

d 2 q (g + r)_ 1 1 
(2tt) 2 (q -r)_ |g s | a ~~ ~tt(2 - a)|r s | a " 2 ' 



(D.7) 



(D.8) 



(D.9) 



where a is suitably greater than a positive number to cut the UV divergence. This can 
be obtained in a similar manner to the formulas (2.17-19) derived in Hl2| . Employing the 
formula twice we find 

d 2 q' (p + q% 1 



A' 



(2ny(p-q>)„\qi\^ 
A 2 



n(l-e)\q 



I I e— 1 



d 2 q' (p + q')- 1 

n{l-e)JJ (2tt) 2 (p -<?')- K 
A 2 \ ( 1 



2e 



2 



X 



7r(2-2e)\ Ps 



2e-2 



P 



87T 2 



2 +0(e). 
r 

Similarly we compute the third integral, 4tt 2 \ 2 J 2 , where 
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d 2 q (p + q)- f d^~ e) q 3 1 



.4 



A x 
2vr 



(2n) 2 (p-q)-J (2n)^q 2 + q 2 
d 2 q (p + q)- 1 



{2izY{p-q)_\q 
1 

tt(1 -e)^^" 1 
0(e). 



|l+e 



(D.ll) 



Since 2J 2 + J 2 = with e — > by analytic continuation, the second and third terms in 
( P-3| ) cancel. 

Finally we turn to the fourth term. By using the solution ( 2.65|) , the fourth term 
becomes — 2NX i f J 4: , where 

WlViTT (D12) 

J (2n) 3 q 2 s +q 2 

Following the prescription in Appendix |C| we regularize this integral as follows. 



3a 



d 2 q fd^q 3 ^qJT 



(2tt) 2 J (27r) 1 ^g 2 + m 2 + g 2 ' 
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(D.13) 



where we insert a virtual mass in 2-plane to cure IR divergence. Performing integration 
we find 



J d = A 



d 2 q 



( 2?r ) 2 yjq 2 a + m 2 



l+e 



-A 



m 



2-e 



2tt(2 



(D.14) 



which vanishes under m — > 0, e — > 0. 

As a result we complete checking the equation (p.3[). 



D.3 Integration formulas 

In this appendix we collect detailed calculation used in this paper. 
We first show the following integration formula. 



sE/ (^( W T^TlP = -^ log 2smh( 2 > ' (D15) 



We regularize the left-hand side such that 



( ..h., m(E Ti)4^r^/^ 71 — i-— . (D.i6) 



P V A. P*0 V (2tt)-< (2f ) 2 + ? + ,J + A/ 2 

By using the formula 



ly ! = £25<£) (D17) 



where h is some positive constant, we can perform the summation. 

(i.h. s . in q * dq * r d q coth ^ — • 



(2tt) J (2vr)- e 2v/g 2 + q 2 s + M 2 



, „ ,_ e * / + i //V<? 2 +<?s 2 + M2 \ 1 i 
g s ctg s I d q I coth(— ^ — ^ ) — 1 1 



(2tt) J (27r)- e V 2v/g 2 + g 2 + M 2 2^<f + g 2 + M 2 , 

(D.18) 

In the second equality, we separate the divergent part as the second term so that the 
first term is convergent. Therefore the first term can be evaluated without calculating the 
integral of extra dimension and results in 



X (2ir) 2^9 2 + M 2 4ir/3 
The second term is calculated as 

f°° q s dq s f d-*g 1 = 1 [°° q s dq s 1 £(|+f) 2 2 

p (2tt) 7 (27r)- 2v /g 2 + g 2 + M 2 2 i p (2vr) (4vr)¥ r(|) Ws J 

= 1 1 r(| + |) / (p 2 + M 2 )^ 
"2( 47r )¥ r(|) I 27r(l-e) 



(D.19) 



I_ e 
' 2 2 



(D.20) 
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Note that e has to be greater than one to make this integral convergent. By analytic 
continuation from some positive number e to zero, we finish evaluating the second term as 

— ^ p ^ M - . Summing up these we obtain (|D.15| ). 
By setting p = for (D.15) we obtain 



R L^i 



d 2 q 



(2tt) 2 (Iff + q 2 + M 2 



2vr/3 



log I 2 sinh( 



(3\M\ 



(D.21) 



We can obtain a similar formula to ( p.21| ) by changing it from integer sum to half- 
integer sum. 



-y 



d 2 q 



1 



(2vr)2 ( 2^4) )2 + gs2 + M2 



27T/3 



log ( 2 cosh( 



(3\M\ 



(D.22) 



One can show this formula by repeating the same calculation by exchanging sinh X, coth X 
into coshX, tanhX, respectively. This is because the summation formula (|D.17f) is corre- 
spondingly changed by 

1 E^^ = ^. (-) 







neZ 



2h 



Next we generalize the formula ( P-22|) in the following manner. 



R ^ 







d 2 q 



log ( 2cosh( v 2 
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(log(2cosh(^))) 
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(D.24) 



The case with a = is obtained above. The procedure to calculate is the same above. We 
regularize the left-hand side such that 



1 



(l.h.s. in (EH = - 



d e q 



log ( 2cosh( v 2 



d 2 q 

J2^J (27P - (2E2)2 + q 2 + ql + M 2 
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Employing the formula ( |D . 1 7| ) we can carry out the sum 

log 2cosh( v 2 J 



(l.h.s. in (OH 



d q f d e q 



coth ( — Y 2 ) 



(2vr) 2 J (2tt)-' 
d 2 q f d~ e q 
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log(2cosh(^|±^; 



coth(— Y ^ 



2 



+ 



(2tt) 2 J (2vr)- e 

c/ 2 g /" c/ _e g 



2 v /| 2 Tg 2 TM 2 
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(27r) 2 i (27r)-2v/g 2 + g 2 + M 2 ' 
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Figure 4: The bosonic 1PI self energy diagram is drawn. The dotted, undotted and wavy 
lines respectively represent scalar, fermion and gauge field. The black and gray bubble 
respectively stand for the planar propagators of scalar and fermion. 



In the second equality we separate the divergent part. Now the first term is convergent, so 
it can be evaluated without taking into account the integral of extra dimension as follows. 



d 2 q 



log ( 2cc.sh(^p^V coth(^p^) - ( ^f±^ 



(2tt) 2 2y/q 2 a + M 2 

-(log(2cosh(^))) a+1 1 //JVM^ 1 



27r/3(a+ 1) 4tt \ 2 J a + 1 

The second term is calculated as 
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(D.28) 

Note that e has to be greater than a + 1 to make this integral convergent. By analytic 
continuation from some positive number e to zero, we finish evaluating the second term as 
— 4^ (f ) M a+i • Summing up these we obtain ( |D.24j) . 



E Diagrammatic analysis 



In this appendix we check the gap equations fl2.58|) and the exact effective action fl2.59|) 
from a diagrammatic point of view. 

First we consider the one-particle irreducible (1PI) diagram for bosonic self energy.^ 
The bosonic 1PI self energy diagram is drawn in Fig.|]. Note that gauge self interaction is 
not involved in these diagrams due to our gauge choice. This diagram encodes an equation 



15 Remind that a 1PI diagram is such that it cannot be split into two by cutting any line from it. 
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Figure 5: The fermionic 1PI self energy diagram is drawn. 



such thatQ 

d?q 



- Nt Sw?Sf {p + q ' r Wrkw G ^ p - q ' ] ¥Tk¥) G%(p - q){r + " r 

d 3 q d 3 q' 1 1 



3A 6 iV 



(27r)3(27r)3g2 + Si3 ( g ) q v + Z B (q>) 



f-tr ,..,„. ' 73 



(2tt) 3 V z7^ + Sf(?) 
where G^ u is the relevant coefficient of the gauge propagator^ given by 

G^ip) = — (^ + <f„, 3 - ^ >3 ^+). (E.2) 

One can easily reproduce the gap equation of boson in fl2.58|) from this equation^. 

Next the fermionic 1PI self energy diagram is shown in Fig|| We can read off an 
equation from this diagram such that 

-E F (p) = N [ ^-irfG^ip -q) * i>f - A 4 iV / f q I (E.3) 

It turns out that this agrees with the gap equation of fermion in ( |2.58| ). 

Finally we check that the effective action (|2.59|) encodes the connected vacuum graph 
in a perturbative way. For this purpose, we first expand —S e ff ( |2.59j ) by — — Sp, which 
is diagrammatically drawn in Fig|| Plugging ([E.lfj and ( E.3|) into the equation expanded 
by — £5,— S^, we obtain the perturbative expansion in terms of X,X4,xq. Each term 
represents a planar diagram contributing to the vacuum graph. These planar diagrams 
are categorized into three classes by matter content: One includes only scalar, another 
only fermion, the other both. The first class, which include only scalar and gauge field 
propagators, is the following.^ 



16 Note that 1PI self energy diagrams are given by — Ss, — £f in our notation. 
17 The gauge propagator is given by (A%(-p')A b v (p)) = S ab {2n) 3 S 3 {p - p l )G IMU {p). 
18 We note that the first term in (E.f ) vanishes. 

19 If we normalize A6 — > then the 6th and 7th terms respectively have the coefficients i and 3. 
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Figure 6: We draw the perturbative expansion of —S e ff in terms of — — ~E F diagram- 
matically, which is expected to describe the connected vacuum graph. 



The 2nd class, which involves only fermion and gauge propagators, is 




The third class which includes both scalar and fermion lines is given by 

- OO n OCQ -i COO + • 



These contributions can be explicitly checked by perturbative calculation of the con- 
nected vacuum graph by using the original action ( [2.71) including gauge field. In fact, one 
can easily see that each diagram has the correct symmetry and sign factors.0 



F Supersymmetric Chern-Simons matter action 

In this appendix we present M = 1,2 Chern-Simons theory coupling to a fundamental 
chiral multiplet in our notation. 



F.l AT = 2 case 

Af = 2 Chern-Simons-matter Lagrangian we use consists of Af = 2 gauge multiplet 



(Ap, a, A, A, .D), and a chiral multiplet (0, if, F). The Lagrangian is |35 



rM=2 _ rN=1 , rH=1 /p i ^ 

*- *-cs i '"matter' K r - 1 J 



The graph including np fermion loops has the sign factor (— l) 1 
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where 



rJV=2 
*~'cs 



%e^{A p d v A p - jA p A u A p ) - AA + 2Da 



•N=2 
matter 



D^D^ + i\){JP + <r)ip + ip\4> + <j>\if) + 4>(a 2 - D)<j> - FF. 



(F.2) 
(F.3) 



Here e 123 = 1,1/) = 7 M -D„. Supersymmetry transformation rule for a vector multipelet is 



A e A p 



A e cr 

A e A Q 
A e \ a 
A e D 

For a chiral multiplet, 



V2 



(e7 M A - ej p X), 
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V2e llp \(^ up F^ + [£>> <r]) - V2e a D, 
V2e 7lp t a (- l -e^F^ + [£>», a]) - y/2e a D, 
--L(e 7 ^ 7 /jA^ + e 7 + e[A, a] + e[a, A]). 



A e 
A e ?/; Q 

A e F 

A e 
A e ^ Q 

A f F 



-V2e<ip, 

y/2{e p ^ a (j) + e a a^-e a F), 
-V5(e^ + eo^ + eA0), 



e„0cr — e 



-V2(e W + e^o" + 0eA). 
Integrating out D, a, A, A, F we obtain M = 2 Chern-Simons action 



?N=2 



27 
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2N 



2k 

and supersymmetry transformation rule 
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-y/2ei/j 



(0T a e7^-e 7 ^», 
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A e 



A' 



2k 



(F.4) 

(F.5) 

(F.6) 
(F.7) 
(F.8) 



(F.9) 
(F.10) 
(F.ll) 
(F.12) 
(F.13) 
(F.14) 



(F.15) 

(F.16) 
(F.17) 
(F.18) 
(F.19) 
(F.20) 
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F.2 M = 1 case 



Af = 1 Chern-Simons-matter Lagrangian we use consists of Af = 1 gauge multiplet (A^, A'), 
where A' is a majorana fermion, and matter multiplet (0, if), where is a complex boson. 
The Lagrangian with a general superpotential W(4>4>) is given by [BB[] 



Af=l 



matter ? 



(F.21) 



where 



r N=l 



matter 



is^»{Apd v A p - jA^Ap) - X' a X' 



(F.22) 



D^D^ + + #A> - 20AV + (00) IF'(00) 2 - 



"( ^ + (^0)(#)) + (#)(^)V''(#) 



Here IF' (re) = =: Vt^. Supersymmetry transformation rule is 
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(F.24) 
(F.25) 

(F.26) 
(F.27) 
(F.28) 
(F.29) 



Integrating out gaugino A' we obtain Af = 1 Chern-Simons action 
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and supersymmetry transformation rule as 



A t A; 



A4 
A e ^ 

A4 
A £ ^ c 



■(0T a e 7 ^-e 7 ^T^ 



^2k 
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(F.32) 
(F.33) 
(F.34) 
(F.35) 
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Now, let us impose the conformal symmetry on this action, which requires the potential 
W(4>(f)) to be quadratic. 

— W — 

W(00) = - — (00) 2 , (F.36) 
where w is an arbitrary number. Then the above action becomes 
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and supersymmetry transformation rule 
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= ^-(0T a e 7 ^-e 7 ^T>), 
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= ^^^0+^0(00)), 
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= ^(€^^0+^6,(00)0). 
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Note that M — 1 SUSY in this theory is enhanced to M = 2 when 



w 



This 



can be explicitly seen by comparing M = 2 SUSY action constructed in Appendix [P~1 
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